fs = 200;
band = [30 50];
data = rand(1,200);
refe = rand(1,200);

figure;plot(data);
figure;plot(refe);


data1 = data-refe;
fdata1 = lgt_firfilter(data1,fs,band,'hann');
figure;plot(fdata1);


fdata2 = lgt_firfilter(data,fs,band,'hann') - lgt_firfilter(refe,fs,band,'hann');
figure;plot(fdata2);